clear;clc;
n=0;
%tic;
%ticBytes(gcp);
m=10250000
parfor i=1:m
    s=randi([0,102500],2,1);
    if(gcd(s(1),s(2))==1)
        n=n+1;
    end 
end
%tocBytes(gcp)
p=n/m;
(6/p)^(1/2)
%toc